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Abstract — We describe how a virtual node abstraction layer 
can be used to coordinate the motion of real mobile nodes 
in a region of 2-space. In particular, we consider how nodes 
in a mobile ad hoc network can arrange themselves along a 
predetermined curve in the plane, and can maintain themselves 
in such a configuration in the presence of changes in the 
underlying mobile ad hoc network, specifically, when nodes may 
join or leave the system or may fail. Our strategy is to allow the 
mobile nodes to implement a virtual layer consisting of mobile 
client nodes, stationary Virtual Nodes (VNs) for predetermined 
zones in the plane, and local broadcast communication. The 
VNs coordinate among themselves to distribute the client nodes 
between zones based on the length of the curve through those 
zones, while each VN directs its zone's local client nodes to move 
themselves to equally spaced locations on the local portion of 
the target curve. 

Index Terms — Motion coordination, virtual nodes, hybrid 
systems, hybrid I/O automata. 

I. Introduction 

Motion coordination is the general problem of achiev- 
ing some global spatial pattern of movement in a set of 
autonomous agents. An important motivation for studying 
distributed motion coordination, that is, coordination among 
agents with only local communication ability and therefore 
limited knowledge about the state of the entire system, 
stems from the developments in the field of mobile sensor 
networks. Previous work in this area includes different co- 
ordination goals, for example: flocking [9], rendezvous [1], 
[10], [13], deployment [2], pattern formation [15], and ag- 
gregation [7]. Owing to the intrinsic decentralized nature 
of sensor network applications like surveillance, search and 
rescue, monitoring, and exploration, centralized or leader 
based approaches are ruled out. However, the lack of central 
control makes the programming task quite difficult. 

In prior work [3], [5], [6], [4], we have developed a notion 
of "virtual nodes" for mobile ad hoc networks. A virtual 
node is an abstract, relatively well-behaved active node that 
is implemented using less well-behaved real nodes. Virtual 
nodes can be used to solve problems such as providing 
atomic memory [5], geographic routing [3], and point-to- 
point routing [4]. 

In this paper, we explore a framework for using virtual 
nodes to solve motion coordination problems. We consider 
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virtual nodes associated with predetermined, well-distributed 
locations in the plane, communicating among themselves and 
with mobile "client nodes" using local broadcast. Roughly 
speaking, each virtual node is designated a certain zone in the 
plane — disjoint from the zones of the other virtual nodes — 
and is emulated by all the mobile nodes that are present 
in its zone. The VN abstraction makes programming easier 
by providing a centralized controller with reliable storage (a 
virtual node) for each disjoint zone in the plane. 

In this paper we describe a framework for using vir- 
tual nodes to solve a simple motion coordination prob- 
lem, namely, uniformly positioning the mobile nodes on 
a known differentiable curve. The framework can be used 
to implement distributed algorithms for more complicated 
motion coordination problems as well. Using the virtual 
node abstraction for solving a coordination problem reduces 
to achieving coordination among stationary centralized con- 
trollers with a priori known locations. Further, writing the 
coordination algorithm for the virtual nodes amounts to 
"plugging in" control functions in the virtual node programs. 

In addition to describing the framework for using virtual 
nodes to solve coordination problems, we alsobriefly describe 
one way of implementing virtual nodes using the real mobile 
nodes. We use the Hybrid I/O Automata (HIOA) mathemat- 
ical framework [12] for describing the components in our 
systems. 

The paper is organized as follows: Section II describes the 
underlying mobile network. Section III describes our virtual 
node layer. Section IV defines the motion coordination 
problem we consider. Section V describes an algorithm for 
solving this motion coordination problem using the virtual 
node layer. Section VI gives the proofs of correctness of the 
algorithm. Section VII presents simulation results for our 
algorithm. Section VIII outlines one way to implement the 
virtual node layer, and Section IX concludes. 

II. The Physical Layer 

Our physical model of the system consists of a finite 
but unknown number of communicating physical nodes in 
a bounded square B in R 2 . We assume that each node has a 
unique identifier from a set I. Formally, our physical layer 
model consists of three types of HIOA (see Figure 1): (1) 
automata PN \ to model physical nodes with identifiers i <G 
X, (2) a LBcast automaton that models the local broadcast 
communication service between the physical nodes, and (3) a 



"real world" automaton RW to model the physical locations 
of all the nodes and the real time. 



Xj,i € X, realtime 




Fig. 1. The Physical Layer: PN automata communicate with each other 
through an LBcast service and receive time and location information 
continuously from RW. 



Figure 2 shows the required components of each automa- 
ton PN \\ it may have other internal variables (initially set 
to unique initial values) and actions, which are not specified 
here. PNi continuously receives from RW the current time 
as the input variable realtime and its position as the input 
variable x^, and communicates its velocity to RW through 
the output variable v^. The speed of PNi is bounded by v c . 
The trajectories of the continuous variable v; and the effects 
of the send and receive actions are unspecified. At each 
point PNi is either in active or inactive mode; we assume 
that, initially, finitely many nodes are active. The fail^ input 
action sets the mode to inactive and all internal variables 
to their initial values, and the recover^ input action sets the 
mode to active. In inactive mode, all internal and output 
actions are disabled, no input action except recover^ affects 
the internal or output variables, and during trajectories, the 
locally-controlled variables remain constant and the velocity 
Vj remains zero. Thus, we assume that, in inactive mode, 
PN i stops moving. We model the departure of a node from 
B as a failure. For convenience, we assume that transitions 
are instantaneous. 

The PNs communicate using a local broadcast service, 
LBcast, which is a generic local broadcast service param- 
eterized by a radius R p and a maximum message delay 
dp. The LBcast(R p , d p ) service guarantees that when PNi 
performs a send(m); action at some time t, the message 



is delivered within the interval [t, t + d p ], by a receive {m)j 
action, to every PN j that remains in active mode and within 
R p distance of PNi for the entire interval [t, t + d p ]. 

The RW automaton (see Figure 3) erves to model realtime 
and locations of all the nodes. It reads the velocity output v, 
from each PNi, i e X, and computes the position Xj from 
the velocity v, for PN \ and the LBcast automaton. LBcast 
requires the node position information because it guarantees 
delivery only between "nearby" nodes. RW also produces 
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Input 




Input fail; 


receive(m)i 




Effect 


fail. 




v,^0 


recover; 




mode <— inactive 


Output 




Other internal variables <— initial 


send(m)i 




Input recover^ 


Variables: 
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mode <— active 
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realtime £ R— ° 
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Internal 






mode € {active, 


inactive} 




Finite set of other variables , 


initially set to unique initial values. 



Fig. 2. Hybrid I/O Automaton PN t . 



Variables: 
Input 

v; 6 R 2 , for each i £ X 
Output 

x; 6 B, for each i€l 

realtime £ R—° 
Internal 

x; £ B, for each i £ X, initially arbitrary 

clock £ R-°, initially 

Trajectories: 
Invariant 

x; £ B, for each i £ X 
Evolve 

x; = Xj, for each i £ X 
d(5Li) = v;, for each i £ X 
realtime = clock 
d(clock) = 1 



Fig. 3. RW automaton. 



realtime for all physical layer components. 

III. The Virtual Node Layer 

The bounded square B is partitioned into a finite set of 
zones Bh, h e H. For simplicity we assume B is a m x m 
square grid, with each grid square corresponding to a zone 
and having sides of length b. Each boundary point of a square 
is unambiguously assigned to one zone. The index set Ti is 
the set of coordinates of the centers of all squares. For each 
Bh, the set Nbrsh contains the zone identifiers of the north, 
south, east, and west neighboring grid squares. 

Our virtual layer abstraction (see Figure 4) consists of: 
(1) client node automata CN \ with identifiers i G X, (2) 
one stationary virtual node automaton VN h for each h € 
TC, located at the center o^ of the square Bh, (3) a virtual 
communication service, VLBcast = LBcast (R v , d v ), for the 
VTVs and the CNs, and (4) an automaton RW to model the 
physical locations of all the CNs and the real time. 

A client node automaton CN i, i € X, is a portion of a 
PN i automaton that has the input variables realtime and 
Xj from the RW automaton and an output variable Vj to the 



.i.i £ X, realtime , 



RW (real world) 




receive(m). 



Fig. 4. Virtual Node Layer: VWs and CNs communicate using the 
VLBcast service. 



RW automaton. With respect to failures, an automaton CTVj 
behaves the same as PN i- CN \ also has send and receive 
actions for interacting with the VLBcast service. 

A virtual node automaton V7V \, h E H, is an MMT 
automaton [11], [14] parameterized by a time upper bound, 
Ammt', it has no realtime clock variable. MMT automata are 
discrete I/O automata that have a "task" structure, which is an 
equivalence relation on the set of locally-controlled actions, 
such that from a point in an execution where a task becomes 
enabled, within at most time duMT, some action in that 
task must occur. VNh can experience a fail/j input, disabling 
internal and output actions, preventing any inputs other than 
recover^ from resulting in state changes, and setting the 
automaton to an initial state. If a recover^ occurs at a failed 
VN, the VN actions become enabled with all tasks restarted. 
If VN h is failed and a CN is in Bh and remains active in 
the zone for d r time, then a recover^ occurs within that d r 
time. V/V/j communicates with other VTVs and CNs using 
the VLBcast service through send/j and receive^ actions. 

VLBcast is an LBcast service (as described in the phys- 
ical layer) for the virtual layer, parameterized by radius R v 
and maximum message delay d v , where R v > b. It allows 
VN \ to communicate with each VN g such that g E Nbrsh, 
and with CNs that are located in Bh- It does not allow CN 
automata to communicate with one another. 

The RW automaton in the virtual layer is similar to the 
one in the physical layer, but here it communicates (through 
the realtime and x variables) only with the CN automata 
and the VLBcast automaton, and not the VN automata. 

This virtual layer will be used in Section V to implement 
a solution to the distributed motion coordination problem. 
Details of how this virtual layer can be implemented using 
the physical layer are in Section VIII. There we further 
discuss the relation between the parameters duMT, d r , d v , 
and R p , the physical layer broadcast radius. 



IV. The Motion Coordination Problem 

A differentiable parameterized curve T is a differentiable 
map P — > B, where the domain set P of parameter values is 
an interval in the real line. The curve T is regular if for every 
p E P, |r'(p)| y^ 0. For a,b E P, the arc length of a regular 
curve r from a to b, is given by s(T,a,b) — J \T'(p)\dp. 
r is said to be parameterized by arc length if for every 
p E P, |r'(p)| = 1. For a curve parameterized by arc length, 
s(r,a, b) = b — a. 

For a given point x € B, if there exists p E P such that 
F(p) = x, then we say that the point x is on the curve F; 
abusing the notation, we write this as x E T. We say that F is 
a simple curve provided for every x E T, r _1 (x) is unique. 
A sequence xi, . . . ,x„ of points in B are said to be evenly 
spaced on a curve T if there exists a sequence of parameter 
values p\ < p2 ■ ■ ■ < p n , such that for each i, 1 < i < n, 
T(pi) = x^ and for each i, 1 < i < n, Pi—pi-i = Pi+i—pi- 

In this paper we fix T to be a simple, differentiable curve 
that is parameterized by arc length. Let Ph — {p E P : 
r(p) E Bh} be the domain of T in zone Bh C B. The local 
part of the curve T in zone Bh is the restriction Th : Ph — > 
Bh- We assume that Ph is convex for every zone Bh C B; it 
may be empty for some Bh- We write \Ph\ for the length of 
the curve Th- We define the quantization of a real number 
x with quantization constant a > as q a {x) = \—~]cr. For 
the remainder of the paper we fix a and write qh as an 
abbreviation for q a (\Ph\)- We write q m i n for the minimum 
nonzero qh, and q max for the maximum qh- 

Our goal is to design an algorithm that runs on the physical 
nodes such that, if there are no failures or recoveries of 
physical nodes after a certain point in time, then: (1) within 
finite time the set of nodes in each zone Bh, h ETC, becomes 
fixed, and the size of the set is "approximately" proportional 
to the quantized length qh, (2) within finite time all physical 
nodes in Bh for which qh ^ are located on Th, and (3) in 
the limit all the nodes in each Bh are evenly spaced on Th- 

V. Solution Using Virtual Node Layer 

The Virtual Node abstraction is used as a means to 
coordinate the movement of client nodes in a zone. A VN 
controls the motion of the CNs in its zone by setting and 
broadcasting target waypoints for the CNs: VN \, h E H, 
periodically receives information from clients in its zone, 
exchanges information with its neighbors, and sends out a 
message containing a calculated target point for each client 
node "assigned" to zone Bh- Informally, VN \ performs two 
tasks when setting the target points: (1) it re-assigns some of 
the CNs that are assigned to itself to neighboring VNs, and 
(2) it sends a target position on T to each CN that is assigned 
to itself. The objective of (1) is to prevent neighboring VTVs 
from getting depleted of CNs and to achieve a distribution of 
CNs over the zones that is proportional to the length of F in 
each zone. The objective of (2) is to space the nodes evenly 
on T within each zone. A CN, in turn, receives its current 
position information from RW and its target location from a 
VN, and continuously computes a velocity vector that will 
take it to its latest received target point. 



In our algorithm each virtual node VN \ uses only in- 
formation about the portions of the target curve T in zone 
Bh and the neighboring zones. For convenience, we assume 
that all client nodes know the complete curve T; we could 
instead model the client nodes in B^ as receiving inputs from 
another automaton about the nature of the curve in zone B^ 
and neighboring zones only. 

A. Client Node Algorithm 

The algorithm for the client node CN(8)i, i El, appears 
in Figure 5. The client node follows a round structure, where 
rounds begin at times that are multiples of S. Recall that 
V7V automata do not have access to realtime whereas CN 
automata do. To help VNs follow the round structure, the 
CNs send "trigger" messages to prompt the VTVs to perform 
transitions. 

At the beginning of each round, a CN sends a cn-update 
message to its local VN (that is, the VN in whose zone 
the CN currently resides). The cn-update message tells the 
local VN the CN's id, its current location in B, and current 
round number. 

The CN then sends an exchange-trigger message d v + e 
later to its local VN, An additional duMT + 2d„ + e time 
later, the CN sends a target-trigger message to its local 
VN. Both these messages are trigger messages that include 
the CN's current location and the current round number, 
used by the local VN to determine whether the CN is in its 
zone and what the current round number is. 

CN i processes only one kind of message, target-update 
messages sent by its assigned VN. Each such message 
describes the new target location x* for CN \, and possibly an 
assignment to a different VN. CN \ continuously computes 
its velocity vector Vj, based on its current position X; and its 
target position x*, as v, = f c (xj — x*)/||xj — x*||, moving 
it with maximum velocity towards the target. 

B. Round structure 

The VN h, h E Ti, algorithm follows the CNs' round 
structure. However, VNs do not have access to the realtime 
variable and must instead rely on trigger messages from 
CNs to determine when enough time has elapsed to perform 
required actions. Here we explain how we implement the 
round structure for a VN. 

Recall that at the beginning of a round, each CN sends 
a cn-update message to its local VN. The CNs then send 
exchange-trigger messages d v + e after the beginning of 
the round, enough time that the cn-update messages have 
already been delivered, signaling to the VN that it has 
received all cn-update messages that were transmitted at 
the beginning of the round in its zone. The VN waits before 
using information from the cn-update messages until it 
receives one of the CNs' exchange-trigger messages. The 
VN then sends vn-update messages to its neighbors. 

Each CN sends a target-trigger message to its local 
VN an additional duMT + 2d v + e time after it sends an 
exchange-trigger message. This additional time is enough 
for all the following to have happened: (1) each neighboring 



Signature: 
Input 2 

receive(m)i, m e ({target-update} x B) 

Output 4 

send(m)i, m e ({cn-update} x X x B x N) 

U ({exchange-trigger, target-trigger} x B x N) <-, 
Internal 
initi 8 

Variables: 10 

Input 

Xj S B 12 

realtime £ R—° 

Output 14 

Vj 6 R' 2 , velocity vector 
Internal i6 

x* £6U{1}, target point, initially _L 

round, next-exch, next-target GNU {-L}, initially _L is 

Transitions: 20 

Internal initi 
Precondition 22 

round = _L 
Effect 24 

round, next-exch, next-target <— [realtime/&~\ 



Input receive ((target-update, target)) t m 

Effect 

if target(i) 7^ null then 30 

x* 4— target(i) 

32 

Output send( (cn-update, i, x,, round))i 

Precondition 34 

realtime = round ■ 5 
Effect 36 

round «— round + 1 

38 

Output send( (exchange-trigger, x;, next-exch))i 

Precondition 40 

realtime = next-exch ■ 5 + d v + e 
Effect 42 

next-exch «— next-exch + 1 

44 
Output send( (target -trigger, xj, next-target)) , 
Precondition a 

realtime = next-target ■ 8 + dj^MT + 3d v + 2e 
Effect 48 

next-target <— next-target + 1 

50 
Trajectories: 
Evolve 52 

if (x; = x* or x* = _L) then v; = 

else Vj = v c ■ (x* — x;)/| |x* — x| 54 

Stop when 

round = _L or realtime = round ■ 5 56 

or next-exch-8 + d v + e or next-target-8 + dMMT + 3d« + 2e 



Fig. 5. Client node CN(8)i automaton. 



VN has received an exchange-trigger message from a 
CN in its zone (d v time), (2) each neighboring VN has 
performed a vn-update transmission to its neighboring VNs, 
including this one {cImmt time), and (3) the neighboring 
VN vn-update messages have arrived (d v time). When a 
VN first receives a target-trigger message for a particular 
round from any CN in its region, it knows it has received any 
vn-update messages from neighboring VNs for the round. 
The VN then performs some computation and transmits a 
target-update message to CNs local to it. 

A target-update message might not be received by a CN 
until d,MMT + 2d v time after the CN sent the target-trigger 
message. This accounts for: (1) the time it can take for the 
target-trigger message to be received by the VN (d v ), (2) 
the time it can take for the VN to perform the target-update 
broadcast (dMivir), and (3) the time for the broadcast to 
be delivered at the CN (d v ). Given the maximum distance 
between a point in one zone and the center of a neighboring 
zone, V2.56 = \J(2>b/2) 2 + (fe/2) 2 , and a constant speed of 
v c for each client node, it can take up to v 2 - 5b time for the 
CN to reach its target. Also, after the CN just arrives in the 
zone it was assigned to, up to y/lOb/3 = V2.5& • | distance 
from where it started, it could find that the local VN is failed, 
in which case it could take up to the d r VW- startup time for 
the VN to recover. 

To ensure a round is long enough for a client node to 
send the cn-update, exchange-trigger, and target-trigger 
messages, receive a target-update message, arrive at its new 
assigned target location, and be sure a virtual node is alive in 
its zone before a new round begins, we require that 6 satisfy 
S > IduMT + 5d v + 2e + max(\ / 2.5b/v c , \/l0b/3v c + d r ). 

C. VN algorithm 

The algorithm for virtual node VN(e,p\,p2)h< h £ H, 
appears in Figure 6, where e <G Z + and pi,p2 € (0,1) 
are parameters of the automaton. VN '/, collects cn-update 
messages sent at the beginning of the round from CNs 
located in its zone, aggregating the location and round 
information from the message in a table, M. When VN \ 
first receives an exchange-trigger message for a particular 
round from any CN in its zone, VN \ tallies and computes 
from its table M the number of client nodes assigned to it 
that it has heard from in the round, and sends this information 
in a vn-update message to all of its neighbors. 

When VHh receives a vn-update message from a neigh- 
boring VN, it stores the CN population and round number 
information from the message in a table, V. When VN h 
first receives a target-trigger message for a particular round 
from any CN in its region, VN n uses the information in its 
tables M and V about the number of CNs in its zone and its 
neighbors' zones to calculate how many of the CNs assigned 
to itself should be reassigned and to which neighboring VNs. 
This is done through the assign function (see Figure 7) 
which calculates a partial function assign mapping CN 
identifiers to zones that they are assigned to. If the number 
of CNs y(h) assigned to VN h exceeds the minimum critical 



Signature: 
Input 

receive(m) h , m g ({cn-update} x 1 x B x N) 

U ({exchange-trigger, target-trigger} x & x N) 
U ({vn-update} x H x N x N) 
Output 

send (m) h 

Constants: 

In = {g G Nbrs: q g + 0} 

State variables: 

M:I-+Ex N, partial map from CN ids to current location and 

round number, initially 0. Accessors: loc, round. 
V:H^NxN, partial map from VN ids to the number of CNs, and 

round number, initially {(g, (0, 0)}} for each g G Nbrs U {h}. 

Accessors: num, round, 
send-buffer, queue of messages, initially 0. 
vn-done, target-done £ Z, initially 0. 

Derived variables: 

locM = X(i G id(M)). loc(M(i)) 
y = \(g G Nbrs U {h}). num(V(g)) 

Transitions: 
Input receive((cn-update, id, loc, wund)) h 
Effect 

if loc G -B/i then 

M^MU {(id, (loc, round))} 

Input receive((exchange-trigger, loc, round)) h 
Effect 

if (loc G -Bh A vn-done ^ round) then 
for each i G id(M) 

if round(M(i)) ^ round then 
M <- M\ {(i,M(i))} 
send-buffer «— send-buffer U {(vn-update, h, \M\, round)} 
vn-done *— round 

Input receive ((vn-update, id, n, round)) h 
Effect 

if id G Nbrs then 
V(id) <— (n, round) 

Input receive ((target -trigger, loc, round)) h 
Effect 

if (loc G -Bfe A target-done ^ round) then 
V(h) <- (\M\, round) 
for each g G Nbrs 

if round(V(g)) ^ round then 
V(g) <- (0, 0) 
let target = calctarget(assign(W(M), y), locM) 

send-buffer <— send-buffer U {(target-update, target)} 
target-done <— round 

Output send(m)^ 
Precondition 

send-buffer ^ Am = head(send-buffer) 
Effect 

send-buffer <— tail(send-buffer) 

Tasks and bounds: 

{send(m)^}, bounds [0, dMMT ] 



Fig. 6. VN(e, p\, p2)h IOA signature, variables, transitions, and tasks, 
implementing motion coordination algorithm with parameters: safety e, and 
damping p\,pi. 



Functions: 

function ass\gn(assignedM: 2 Z , y: Nbrs L){h} — > N): X — > H = 
assign: X — > H, initially {((', A)} for each i £ assignedM 
n: N, initially j/(/i) 
ra: N, initially 
if y(h) > e then 
if g/! 7^ then 

let lower = {g G A»: f^2/C0 > 2/(fl)} 
for each g 6 /ower 

ra <- min(Lp 2 • [^f j/(Ji) - j/(s)]/2(|Wer|+l)J, n - e) 
update assign by reassigning ra nodes from h to g 
n <— n — ra 
else if /n = then 

let tower = {</ £ M>r.s : j/(/i) > y(g)} 
for each g £ /ower 

ra <- min(Lp2 ■ [y{h) - y(g)]/2(\lower\+l)\, n-e) 
update assign by reassigning ra nodes from h to g 
n <— n — ra 
else 
ra< _ L(s/(h) - e)/|M|J 
for each g 6 /n 

update assign by reassigning ra nodes from h to g 
return assign 

function calctarget(asv>7g«: J -> W, tocM: J -> B): J -> B = 

seq, indexed list of pairs in P X X, initially the list, for each i£l: 
assign(i)= h AlocM(i) € F^, of (p,i) where p= T^ (locM(i)), 
sorted by p, then i 
for each j £ X : assign (i) ^ n«H 
if assignii) = g ^ h then 

locM(i) <— o 9 
else if locM(i) fr h then 

locM(i) <— choose {min xe r h {<iist(x, ZocM(i))}} 
else let p = r^ 1 (/ocM(j)), .se^fc) = (p, i) 
if /t = first(se?) then /ocM(i) <- r\(inf(P h )) 
else if /: = last(.veg) then locM(i) <— I\(sup(Ph)) 
else let «;g(fc - 1) = (pk_i,ifc_i), .vea(fc + 1) = (pfc+l,*fc+i) 
locMii) «- F h (p + Pl • ( Pfc - 1 + Pfc+1 - p)) 
return /ocM 



Fig. 7. VN(e, pi, p2)h IOA functions. 



number e, then the assign function reassigns some of the 
CTVs to neighbors of VN h . 

Let /rift, denote the set of neighboring VNs of VNh that 
are on the curve T and yh(g), 9 € Nbrsh U {ft.}, denote 
the number num(Vh(g)) of CJVs assigned to VW 9 . If gft, 7^ 
0, meaning WV/, is on the curve (lines 7-11), then we let 
lowerh denote the subset of Nbrsh that are on the curve and 
have fewer assigned CNs than VNh has after normalizing 
with — . For each g <G lowerh, VNh reassigns the smaller 
of the following two quantities of CNs to VN g : (1) ra = 
92 ■ [^Vh(h) - y h (g)]/2(\lower h \ + 1), where p 2 < 1 is a 
damping factor, and (2) the remaining number of CNs over 
e still assigned to VNh- 

If qh = 0, meaning VNh is not on the curve, and VNh has 
no neighbors on the curve (lines 13-17), then we let lowerh 
denote the subset of Nbrsh with fewer assigned CNs than 
VNh- For each g e lowerh, VNh reassigns the smaller of 
the following two quantities of CNs: (1) ra = p 2 - [yh(h) — 
yh(g)]/2(\lowerh\ + 1) and (2) the remaining number of 
CNs over e still assigned to VNh- 

VNh is on a boundary if qh = 0, but there is a g £ Nbrsh 



with q g 7^ 0. In this case, yh(h) — e of VNh's CNs are 
assigned equally to neighbors in Irih (lines 19-22). 

The client assignments are then used to calculate new 
target points for local CNs through the calctarget function 
(see Figure 7). This function assigns to every CN i assigned 
to VN h a target point locM h {i) <G B g ,g e Nbrs h U {h}, to 
move to. The target point locMh(i) is computed as follows: 
If CN i is assigned to VN g , g 7^ h, then its target is set 
to the center o g of B g (lines 30-31); if CN i is assigned to 
VN ft but is not located on the curve Th then its target is 
set to the nearest point on the curve, nondeterministically 
choosing one if there are several (lines 32-33); if CN \ is 
either the first or last client node on Th then its target is set 
to the corresponding endpoint of Tft (lines 35-36); if CN , 
is on the curve but is not the first or last client node then 
its target is moved to the mid-point of the locations of the 
preceding and succeeding CNs on the curve (line 38). For 
the last two computations a sequence seq of nodes on the 
curve sorted by curve location is used (line 27). 

VNh finally broadcasts the new target waypoints for the 
round through a target-update message to its CNs. 

VI. Correctness of Algorithm 

We say CNi,i e I, is active in round t if its mode is 
active for the duration of round t. A VNh, h E TL, is active 
in round t if there is some active CN \ with Xj e Bh for the 
duration of rounds t — 1 and t. Thus, none of the VNs is 
active in the starting round. We use the following notation: 
In(t) is the set of ids h € Ti of VNs that are active in round 
t and for which qh 7^ 0. Out(t) is the set of ids h € H of 
VNs that are active in round t and for which qh = 0. C(t) 
is the set of active CNs at round t, and c7j„(£) and C ou t(t) 
are the sets of active CNs located in zones with ids in In(t) 
and Out(t), respectively, at the beginning of round t. 

For any pair of neighboring zones B g and Bh, and for any 
round t, we use y g (h)(t) to refer to the value of y g (h) at 
the point in time in round t when VN g finishes processing 
the first target-trigger message of round t it receives. For 
an Y fi9 *= Nbrsh U {h}, in the absence of failures and 
recoveries of CNs in round t, yf(h)(t) — y g (h)(t); we write 
this simply as yh(t). We present a sequence of lemmas that 
together establish the following theorem: 

Theorem 1: If there are no failures or recoveries of client 
nodes at or after some round io. tnen within a finite number 
of rounds after to '■ 

(1) the set of CNs assigned to each VNh, h G H, becomes 
fixed, and the size of the set is proportional to the quantized 
length Oft within a constant additive term — ^ TO ~ ' , and 

° *" qminp2 

(2) all client nodes in Bh for which qh 7^ are located on 
Tft and evenly spaced on Th in the limit. 

For the rest of this section we fix a particular round number 
to and assume that no failures or recoveries of CNs occurs 
at or after round to- The first lemma states some basic facts 
about the assign function (see Figure 7): 

Lemma 1: In every round t > io : (1) If 2/ft.(t) > e for 
some heH, then y h (t + 1) > e, (2) In{t) C In(t + 1), 



(3) Out{t) C0ut(t+1), (4) C m (t) C C in (t + 1), and (5) 
C out (t + l)CC out (t). 

Proof: We fix round t > to- (1) From line 6 of the 
assign function (Figure 7) it is clear that VN h , ft E ri, 
reassigns some of its CWs in round t only if yhit) > e. 

(2) For any VN h , ft E In(t), if y h {t) < e then VN h 
does assign CWs, and y^(t + 1) = j/ft.(£), otherwise, from 
line 16 of Figure 7 it follows that yh(t + 1) > e. In both 
cases ft, € /n(£ + 1). 

(3) Same as (2). 

(4) Consider CN \, i E Ci n (t), such that CN i is assigned 
to VNh, ft E In(t). From lines 7-11 of Figure 7 we see 
that CN i is assigned to some VN g , g E Irih U {ft}. Since 
In h U {ft} C 7n(t + 1), the result follows. 

(5) As there are no failures and recoveries of CNs, C(t) = 
C(t + 1). By definition, C in {t) U C out (t) = C(t), C ln (t) n 
C out (t) - 0, and C m (t + 1) U C out (t + 1) = C(t + 1), 
Ci n (t + 1) n C ou t(t + 1) = 0. The result follows from part 
(4). ■ 

The next lemma states a key property of the assign func- 
tion after round t : VN g , g E Out(t), is never assigned a 
larger number of CNs in round t+1 than the largest number 
of CNs that were assigned to any of VN g 's neighbors in 
round t. Similarly, VN g , g E In(t), never gets a density 



y g (t+i) 

la 



of CNs CNs in round t+1 that is greater than the 



highest density of its neighbors in round t. 

Lemma 2: In every round t > to, for g,h E H and ft E 

Nbrs g : (1) If g, ft E Out(t), y h (t) = max feNbrS!i yf(t), 
and y g (t) < y h (t), then y g (t + 1) < y h (t) - 1, and 



viit) 



Vg(t) 



(2) If g,h € In(t), ^ = maXfeNbrSg y -^, and *£* .' 



Vh(t) 



, then 



2/ g (*+l) < 2/h(*) 



Proof: (1) Fix 5, ft and i, as in the statement of 
the lemma. Since yh(t) > y g (t) and g, ft <G Out(t), 
we see from line 16 of Figure 7 that the number of 
CNs that V7V S is assigned from VNh in round t is at 
most P2(yh(t) — y g (t))/2(\lowerh{t)\ + 1). This is at most 
P2(Vh(t) - j/ 3 (t))/4, because y h (t) > y g {t) implies that 
lower hit) > 1. Then, the total number of CNs assigned 
to VN g in round t by all four of its neighbors is at most 
P2(yh(t)-y g it)). Therefore, y g (t + l)< y g (t) + p 2 {yh{t)- 
y g (tj) = p 2 yh(t) + (1 - P2)y g (t). As pi < 1, we have 
y g {t + 1) < yhif). The result follows from integrality of 
y g (t + l) and y h {t). 

(2) As in part 1, fix g, ft and t. Here ^^ > ^^ and 
g, ft G In(t). From line 10 of Figure 7, it follows that the 
number of CNs that V7V g is assigned from VN \ in round t 
is atmostp2(f£yfc(*)-2/g(*))/2(|Zou;er/ l (i)| + l). This is at 
most p 2 (— yJt)-y g {t))/4:, Then, the total number of CNs 
assigned to VN g in round t by all four of its neighbors is 
at most p 2 i^yh(t) - y g {t)). Therefore, y g (t + 1) < (1 - 



2/ g (*+l) 



li/gffl 



P2)l/ fl (t) + P2^^(*), that is i^^ < (1 - p 2 )^f + 

jrtW A <• x h « g (*+i) < jmW a j j 

' ih ' z ' q g qh F 



calculation shows that if Vh } ' ^ ^-^-, then 



Ha Hh 

Vh(t) I Vgjt) t . y h (t) _ y g (t) 



maxk(t) = 



max 
max 



<ih 



>hi 



<li, 



1g 



The next lemma states that there exists a round T out that 
is reached within a finite number of rounds after to, such 
that in every round t > T out , the set of CNs assigned to 
VNh, ft E Outif), does not change. 

Lemma 3: There exists a round T out > to sucn that in 
any round t > T out , the set of CNs assigned to VNh, ft E 
Out(t), is unchanged. 

Proof: First, we show that the number of CNs assigned 
to VNh, ft E Out(t), remains unchanged, that is yh(t + l) = 
yh(t). Let N out be the total number of ft E H such that 
qh = 0. For any k, 1 < k < iV out , we define maxk(t) to 
be the fc t/l largest number of CWs that are assigned to any 
VTV/j, ft € Out{t), at the beginning of round t > to'. 

■{y h (t) :hEOutit)}, iffc = l 

:{l/ h (t) : ft G Outit) A 

y/i(i) < maxk-i(t)}, otherwise. 

Let maxvnskit) be the set of VN ids that have maxk(t) 
CNs assigned to them. If there exists an I, 1 < I < N out , 
such that Vft E Outit) : maxi(t) > yh(t), then for all k, 
I < k < N out , maxkit) = and maxvnsk(t) = 0. 

Let E(t) = (\C ou t(t)\, maxi(t),\maxvnsi(t)\,. . ., 
maxN ou , t it),\maxviiXN out (t)\). Let w be the minimum 
Vh(to) ror an Y h E Out(to), and S = {h E Outifo) '■ 
Vh(to) — w }- Observe that if w < e, then E m i n = 
{w\S\,w, \S\,0, 0. . . ,0, 0) is a minimum value for E(t), 
otherwise E m i n = (ej<S'|,e, j/S'j,0,0. . .,0,0) is a minimum 
value. It suffices to show that for any round t > to, either 
E(t + 1) = E(t), that is, t = T out , or E(t + 1) is less 
than E(t) by some constant amount, meaning there is a 
k, 1 < k < N out , such that for every 1, 1 < I < k, the 
I th component of Eit + 1) is equal to the I th component of 
E(t), and the k th component of E(t + 1) is less than the 
k th component of E(t) by at least 1. 

Consider any round t after to- From Lemma 1 we know 
that |C «t(t + l)| < \C out (t)\. If \C out (t + l)\ < \C out {t)\, 
then the first component of E(t + 1) is less than that of E{t) 
by at least 1. Otherwise, \C out (t + 1)| = |C out (t)|. If for 
every ft E Outit), ra = for all g E lower hit) (see line 16 
of Figure 7), then none of the CNs in C out (t) are reassigned 
in round t + 1, and E(t+1) = Eit). Setting T out = t, we are 
done. Otherwise, there exists a nonempty set of VNs with 
ids in Outit) that reassign some CNs to a neighboring VN. 
We select the nonempty set A of such VNs with the highest 
number of assigned CWs. Let A C maxvnskit), for some 
k, 1 < k < N^. 

For any g E Outit) with y g (t) < maxk(t), the maximum 
value of yhit) for any ft E Nbrs g such that VN g gets some 
CWs from VN \ in round t is at most maxkit). From Part(l) 
of Lemma 2 it follows that y g (t + 1) < maxkit) — 1. 

For any VNh, h E A, since no VN with y > 
maxkit) assigns any CWs to VNh, yhit + 1) = yhit) — 
2~2 g eiower h (t) ^sW* where ra g is the number of CWs VN h 
assigns to its neighbor VN g in round t. We have shown 
above that for any g E Outit), if y g (t) < maxkit) then 
y g if, + 1) < maxkit) — 1. There are two possible cases: (1) 



if maxvnsk(t) = A, then the k max decreases, maxk{t + 
1) < maxk(t) — 1. That is, the (2fc + l) s * component of 
E decreases by at least 1, and (2) if A C maxvnsk(t), 
then maxk(t + 1) = maxk{t) and \maxvnsk(t 4- 1)| = 
|maxwnsfc(i)| — |A|. That is, the (2fc + 2)" d component 
of E decreases by at least 1. This implies that there exists 
T out , such that the number of CNs assigned to each VNh, 
h E Outit), t > T out , remains unchanged. 

Now suppose the set of CNs assigned to V7V \ changes 
in some round t > T out . Since yhif + 1) = jjh(t) for all 
h E Out{t). Summing, \C out (t + 1)| = \C out (t)\ and using 
Lemma 1 we get C ou t(t+l) = C ou t(t). The only way the set 
of CNs assigned to VN h could change, without changing 
yh and the set C ou t, is if there existed a cyclic sequence of 
VNs with ids in Out(t) in which each VN gives up c > 
CNs to its successor VN in the sequence, and receives c 
CNs from its predecessor. However, such a cycle of VNs 
cannot exist because the lower set imposes a strict partial 
ordering on the VNs, ■ 

For the rest of the section we fix T out to be the first round 
after to, at which the property stated by Lemma 3 holds. 
Lemma 3 implies that in every round t > T out , In(t) — 
In{T out ), Out(t) = Out{T out ), C in (t) = C in {T out ), 
and C outit) = C ou t(T out ); we denote these simply as 
In,Out,Cin, and C ou t- The next lemma states a property 
similar to that of Lemma 3 for VN h, h E In, and its proof 
is similar to the proofs of Lemma 3, and uses part (2) of 
Lemma 2. 

Lemma 4: There exists a round T sta b > T out such that in 
every round t > T sta b, the set of CNs assigned to VN \, 
h E In, is unchanged. 

The following lemma bounds the total number of CNs 
located in zones with ids in Out to be 0(to 3 ). 

Lemma 5: In every round t > T outl \C ou t{t)\ — 0(m 3 ). 
Proof: From Lemma 3, the set of CNs assigned to each 
VN h, h E Outit), is unchanged in every round t > T out . 
This implies that in any round t > T ou t, the number of CNs 
assigned by VNh to any of its neighbors is 0. Therefore, 
from line 20 of Figure 7, for any boundary VN g , iy g (t) — 
e)/\In g \ < 1. In g is the (constant) set of h E Nbrs g with 
q h ^ 0. Since \In g \ < 4, y g {t) < 4 + e. From line 16 
of Figure 7, for any non-boundary VN g , g E Outit), that 
is 1-hop away from a boundary VN h , ^2~{t)\+l) < L 
Since \lower g {t)\ < 4, y g (t) < — +4 + e. Inducting on the 
number of hops, the maximum number of CNs assigned to 
a VN g , g E Outit), at I hops from the boundary is at most 

— + e + 4. Since for any I, 1 < I < 2m — 1, there can be at 

pi j — — 

most to VNs at Z-hop distance from the boundary, summing 

10to 2 (2to-1) _ ^^^,3^ 
P2 



gives \C out \ < (e + 4) (2m - l)m + 



0(to 3 ). 



For the rest of the section we fix T s t a b to be the first round 
after T out , at which the property stated by Lemma 4 holds. 
The next lemma states that the number of CNs assigned to 
each VNh, h G In, in the stable assignment after T sta b is 



proportional to qh within a constant additive term. 
Lemma 6: In every round t > T sta b, for g,h e 7n(i): 



Vh(t) y g {t) 



Qh 



<lg 



< 



10(2to- 1) 



qminP2 



Proof: Consider a pair of VNs for neighboring zones 
B g and Bh, g,h <G In. Assume w.l.o.g. j//j(i) > y g (t). From 
line 10 of Figure 7, it follows that p 2 i^Vhit) - y g (t)) < 

2(\lower h it)\+l). Since \lower h (t)\ < 4, |^ - ^| < 
< — — — . By induction on the number of hops from 1 



QgP2 — q m inP2 

to 2m — 1 between any two VNs, the result follows. ■ 

From line 33 of Figure 7, it follows immediately that by 
the beginning of round T sta b + 2, all CNs in Cj„ are located 
on the curve T, This establishes that the VN algorithm 
satisfies our second goal. The next lemma states that the 
locations of the CNs in each zone Bh, h G In, are evenly 
spaced on Th in the limit. 

Lemma 7: Consider a sequence of rounds t\ = 
T s tab, ■■■ ,t n - As n — > oo, the locations of CNs in Bh, 
h <G In, are evenly spaced on Th- 

Proof: From Lemma 4 we know that the set of CNs 
assigned to each VNh, h E In, remains unchanged. Then, 
at the beginning of round ti, every CN assigned to VN \ is 
located in Bh and is on the curve Th. Assume w.l.o.g. that 
VNh is assigned at least two CNs. Then, at the beginning 
of round t%, one CN is positioned at each endpoint of Th, 
namely at r^(w/(P/ l )) and r^(sup(P/ l )). From lines 35- 
36 of Figure 7, we see that the target points for these 
endpoint CNs are not changed in successive rounds. Let 
seq h it 2 ) = (p ,i(o)) , ■ ■ ■ , (Pn+iJ(n+i)), where y h = n + 2, 
Po = infiPh), and p n +i = supiPh). From line 38 of 
Figure 7, for any i, 1 < % < n, the i th element in seqh 
at round tk, fc > 2, is given by: 



Piitk+i) =Pi(tk) + Pi 



Pi-i(tk) +Pi+i(tk) 



-Pi{tk) 



For the endpoints, piitk+i) = Piifk)- Let the i th evenly 
spaced point on the curve Th between the two endpoints be 
Xj. The parameter value pi corresponding to Xj is given by 
Pi = Po + ^rjiPn+i — Po)- In what follows, we show that as 
n —* oo, the pi converge to pi for every i, < i < n+ 1, that 
is, the location of the non-endpoint CNs are evenly spaced 
on Th- The rest of this proof is exactly the same as the proof 
of Theorem 3 in [8] in which the authors prove convergence 
of points on a straight line with even spacing. 

Observe that k = ±(pVi+pVk) = (l-pi)Pi + ^(p»-i + 
pi+i). Define error at step k, k > 2, as e,(fc) = p,(ife) ~ 
pi. Therefore, for each i, 2 < i < n — 1, &iik + 1) = 
Piitk+i) - Pi = (1 - pi)ei(k) + ^-(ej-i(fc) + e i+ i(fc)), 
ei(fc + 1) = (1 - pi)ei(fc) + ^e 2 (fc), and e n (k + 1) = 
(1 — pi)e n (k) + Qe n -iik). The matrix for this can be written 



as: e(k + 1) = Te(k), where T is an n x n matrix: 



I -Pi 
Pi/2 






Pi/2 
1-pi pi/2 









Pi/2 





Pi 
Pi 






Pi/2 

Pi/2 



Using symmetry of T, p\ < 1, and some standard theorems 
from control theory, it follows that the largest eigenvalue of 
T is less than 1. This implies linik^ooT k = 0, which implies 
lim k ^ 00 e(k) = 0. ■ 

VII. Simulation Results 

In this section we briefly describe the performance of the 
algorithm as observed in Matlab simulations. In particular 
we study the stabilization time T stao , that is, the number of 
rounds required to achieve a stable distribution of PNs over 
zones in 7i, for different types of target curves. The inputs 
to the Matlab function that models the assign function of 
our algorithm are: (1) a value of the parameter e of the 
algorithm, and (2) two m x m matrices corresponding to 
the initial and the target distribution of PNs respectively. 
We have performed simulations for different values of m 
and for different target distributions. In all the experiments, 
we fix the value of e to be 5, the number of participating 
PNs to be 10 5 , and the initial locations of all PNs to be the 
lower right corner of B. 

The bar charts in in Figure 8 show the distribution of 
the PNs (dark bars) at different stages in a typical run of 
the algorithm, and the target distribution (light bars). From 
top to bottom, the charts show the distribution of the PNs 
after the first round, at the end of round T out , and the stable 
distribution attained at round T sta b, respectively. 

The plots in Figure 9 show the values of T out and T stau 
for values of m ranging from 10 to 20, and for three different 
curves. CI covers a rectangular region at the lower-right 
corner of B, C2 covers a rectangular region at the top-left 
corner of B, and C3 is an annular ring at the center of B. 
From these plots we observe that if the curve is located far 
from the initial position of the PNs, as in the case of C2, 
then the stabilization time T sta b is dominated by T out . 

VIII. Implementing the Virtual Node Layer 

In addition to client CN i, a physical node PN i, i s I, in 
zone Bh runs a TOBcastih service and a VNEi t h, h G H, 
algorithm (see Figure 10) to help implement each virtual 
node VN h and the VLBcast service of the virtual layer. 

In this section we present a sketch of our implementation 
of the virtual layer by the physical layer. Our implementation 
is an adaptation of techniques from [4] to emulate a virtual 
mobile node. The only substantive changes made in our 
current implementation are: (1) the changing of virtual node 
locations to be stationary, (2) the replacement of a periodic 
location update with a continuous real-time location update, 
and (3) the restart of a virtual node as soon as a physical 
node discovers it is in a failed virtual node's zone. The 
virtual nodes we implement here are also modeled differently 




Fig. 8. Simulation results show the actual (dark) and target (light) 
distribution of PNs over Ti., at the end of various rounds: the first round 
(top), at T ou t (middle), and at T stao (bottom). 
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Fig. 9. Plot of T ou t and T ata i, for three different kinds of curves with 
varying length m of each side of the bounded plane B. CI covers a 
rectangular region at the lower-right corner of B, C2 covers a rectangular 
region at the top-left corner of B, and C3 is an annular ring at the center 
of B. 




Fig. 10. PNi's subautomata: A physical node runs several programs, 
including VNE and TOBcast automata as well as a CN automaton. 



from those in [4], as MMT automata, rather than simple I/O 
automata. 

We use a standard replicated state machine approach to 
implement robust virtual nodes that takes advantage of a 
TOBcast service to ensure that all VNEs in a zone receive 
the same messages in the same order. Using the LBcast 
service of the physical nodes and common knowledge about 
realtime, the totally ordered broadcast service TOBcast 
for a zone can be implemented as follows: At the time of 
sending, a message is tagged with the sender's identifier, 
zone id, a sequence number, and a timestamp, which is the 
current value of realtime. The tags define a total order on 
sent messages, used in delivery. Before delivering a message 



TOBcasti^h waits until d p + e time has elapsed since it 
was sent, ensuring that earlier messages were received. It 
then delivers the messages for the timestamp in order of 
sender id and sequence number. (As a technical detail, we 
actually order all join-req messages (described shortly) after 
any other messages for a particular time.) TOBcastih only 
processes messages tagged for zone Bh- 

Each VNE. i: h independently maintains the state of VN h 
and simulates performing actions of the VN on that state. In 
order to keep the state replication consistent across different 
VNEs running on different physical nodes in the same 
zone, when VNE^h wants to simulate an action of the VN 
(such as that of receiving a message at the VN that was 
actually received by VNE^h), it broadcasts a suggestion to 
perform the action to the other VNEs of the zone using the 
TOBcast service. When an action suggestion is received by 
VNEi : h, it is saved in a pending-action queue. Actions are 
removed from a pending-action queue in order by VNEi t h 
and simulated on VNEu-Js local version of the VN state. 
A completed action is then moved into a completed- action 
queue, referenced by VNEi : h to prevent reprocessing of 
completed actions. 

When a VNE enters a zone, it executes a join protocol to 
get the zone's VN state. The join protocol begins by using 
TOBcast to send a join-req message. Whenever a VNE 
receives its own join-req message, it starts saving messages 
to process in its pending- action queue. If a VNE that has 
already joined receives the join-req, it uses TOBcast to send 
a join-ack containing a copy of its version of the VN state. 
When the joining VNE receives the join-ack, it copies the 
included VN state and starts processing the actions in its 
pending-action queue. If a VNE's join-req is not answered 
in 3d p + 3e time, indicating the VN is failed, the VNE will 
reset the VN e time later by using TOBcast to send a reset 
message. When a VNE receives a reset message, it sets the 
VN state to its initial state, clears the pending- action queue, 
and starts simulating the VN, 

Theorem 2: Assuming R p > \/5b, the TOBcastih, 
VNEi t h, i € I, h € H, and trivial client implementation 
correctly implement the Virtual Node abstraction with VN 
task upper time bound duMT = 2d p + 2e, VJV-startup time 
d r = Ad p + 5e, VLBcast broadcast radius R v > b, and 
VLBcast maximum message delay d v = 2d p + e. 

Proof: The correctness of the implementation of the 
Virtual Node layer largely follows from the proof of cor- 
rectness for the implementation of the VMN layer in [4]. 
We here discuss the correctness of the implementation with 
respect to: (1) the task upper bound, (2) the VTV-startup time, 
and (3) the requirements for LBcast and VLBcast. 

(1) Once one of an abstract VNh's output or internal tran- 
sitions is enabled, the precondition for sending a suggestion 
to simulate the action through TOBcast is satisfied at all 
VNEi t h for PNi in Bh, and the broadcast occurs. It takes 
at most dp + e time for the message to be delivered at other 
VNEi : h for PNi in Bh, after which the action is simulated. 
However, it is possible for all active VNEs to fail right after 



sending a join-ack to a new VNE and before proposing 
an enabled action, leaving the new VNE to broadcast the 
simulation proposal d p +e later, when it receives the join-ack. 
Given that PN transitions are assumed to be instantaneous, 
diAMT = 2d p + 2e. 

(2) If PN \ enters a zone B^ with a failed VN, its 
VNE^h's join-req will not be answered in 3d p + 3e time, 
and the VNE will send a reset message an additional e later. 
It takes the VNE at most d p + e time to receive the reset 
message and restart the VN. The total time 4d p + 5e for a 
joining node to succeed in restarting a VN is d r . 

(3) As in [4], d v = 2d p + e since the underlying LBcast 
service used to implement VLBcast takes up to d p time 
to deliver a transmitted message from a VN or CN , after 
which TO Beast takes an additional d p + e time to redeliver a 
message at a receiving VN. Also similarly to [4], we require 
that R p > \fbb, in order to guarantee that R v > b, allowing 
a CN i in Bh, i E X, h E H, and VN h to communicate, and 
a VN h (located at Oh) and each of its neighboring zones' 
VN g , g E Nbrs(h), (located at o g ) to communicate. This is 

because a VNE emulating a zone Bh can be as far away as 
y / (2b) 2 + b 2 from a VNE emulating the VN of neighboring 
zone B g . To guarantee the two can communicate while 
emulating their respective VTVs, the broadcast radius R p of 
the physical LBcast service must be be at least y/Sb. Unlike 
[4], however, we do not require an additional tolerance factor 
to account for periodic location updates from the RW; here, 
the RW automaton is assumed to continually update the 
VNE of its current location. ■ 

IX. Future work and extensions 

We believe the framework introduced in this paper can 
be useful in simplifying the coordination of mobile nodes 
to solve a variety of other, more complex, problems. One 
promising avenue of future work would be to employ our 
framework for some of those problems. As one example, in 
the control algorithm presented in this paper, each virtual 
node VN uses only local information about the target curve 
r. This use of only local information should adapt well to a 
problem extension where the curve is dynamically changing. 
The curve (or point, even) could be moving targets being 
tracked. In this case, the framework for coordination of nodes 
we present here is useful for two reasons: (1) maintaining 
alive VTVs to detect targets and (2) guiding physical nodes 
to the moving targets. 

Also of interest is an in-depth analysis of the stability 
of the control algorithm employed in this paper in the 
face of some regular rate of physical node addition and 
removal. While there has been a paucity of work in this 
area, such analyses are extremely important in the evaluation 
of solutions to most any control algorithms for dynamic 
systems. 
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Appendix 



Signature: 
Input 

receive(m); f,, m a client message 

TOBcast-rcv(m)i h> wi a TOBcast message 
Output 

send(m)j ; /j, m a client message 

TOBcast(m)j h m a TOBcast message 
Internal 

zone-update^ 

Joir\ h 

restart; , h 

init-action(art) ih , act £ VN h .sig \ inputs 

simulate-action(art) i h , act e VN h .sig 
ack-join l>fe 

Variables: 
Input 

x; £ B, current location of mobile node 
realtime £ R— ° 
Internal 

status £ {joining, listening, active}, initially active 

h eHU{±}, zone id, initially _L 

val £ VN h .states, state of VW h , initially VN h . start 

pending-join, max id of pending join reqs, initially 

completed -join, max id of answered join reqs, initially 

join-id, time of join-req, initially 

pending -actions, queue of VN^. actions to be simulated, initially 

completed -actions, queue of VN^. actions simulated, initially 

TOBcast-out, queue of outgoing TOBcast msgs, initially 

local-out, queue of outgoing client messages, initially 



Trajectories: 

Stop when any Precondition is satisfied 



Fig. 11. Signature, variables, trajectories of VNEi ^ algorithm imple- 
menting VNh. 



Input receive (m) ith 
Effect 

TOBcast-out <— TOBcast-out U {(simulate, (receive, hi), -L)} 

Output send(m)j ^ 
Precondition 

local-out ^B Am = head(local-out) 
Effect 

local-out <— ta\\(local-out) 

Internal init-action(art) ih 
Precondition 

status = active A x £ B h A 5(va/, art) 7^ _L 
Effect 

TOBcast-out <— TOBcast-out U {(simulate, art, (realtime, i))} 

Internal joinj ^ 
Precondition 

.rtarws = idle Ax £ Bh 
Effect 

status <— joining 

join-id <— realtime 

TOBcast-out <— TOBcast-out U {(join-req, ±, join-id)} 

Internal restart^ 
Precondition 

status = listening A x £ 1?^ A realtime = join-id + 3d p + 4e 
Effect 

TOBcast-out <- TOBcast-out U {(reset)} 

Internal zone-update ijh 
Precondition 

x^Bh 
Effect 

status *— idle 

ft <— id of zone h' such that x £ i?^/ 

vaZ <— VN h . start 

pending -actions <— 

Internal simulate-action(art) i ,/ l 
Precondition 

sfaft« = active Ax £ Bfj Ahead(pending-actions) = (simulate, act, old) 
Effect 

dequeue (pending -actions ) 

if ((simulate, art, oid) £ completed -actions A<5(va/, art) ^ _l_) then 
vaZ <— 5(va/, art) 
if art = (send, m) then 

local-out <— local-out U {m} 
completed -actions <— completed -actions U {(simulate, art, of'd)} 

Internal ack-joinj, fc 
Precondition 

status = active Ai £ i?f, A pending-join > completed -join 
pending -actions = AVart £ VNh.sig \ inputs: 8(val, act)= _L 
Effect 

TOBcast-out <— TOBcast-out U {(J0in-ack,(i<a/, completed -actions) , pending-join)} 
completed -join <— pending -join 

Input TOBcast-rcv((opfype, param, oid))i t f l 
Effect 

if optype = simulate then 
if ,«taft« = listening or active then 

enqueue (pending -actions, (simulate, param, oid)) 
if optype = join-req then 

pending-join *— max(pending-join, oid) 
if (status = joining AcM = join-id) then 
.stents <— listening 
if optype = join-ack then 

completed -join *— max(completed-join, oid) 
if (status = listening and oid > join-id) then 
.vtato <— active 

(va/, completed -actions) <— param 
if optype = reset then 
siato <— active 
pending -actions <— 

Output TOBcast(m) ih 
Precondition 

TOBcast-out ^ Am = head(7"OBca.sf-or<r) 
Effect 

TOBcast-out <- tail(TOBcai(-o«t) 



Fig. 12. Transitions of VNEi^ algorithm. 



